Ring dark solitary waves: experiment versus theory 
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Theoretical and experimental results on optical ring dark solitary waves are presented, empha- 
sizing the interplay between initial dark beam contrast, phase-shift magnitude, background-beam 
intensity and saturation of the nonlinearity are presented. The results are found to confirm quali- 
tatively the existing analytical theory and are in agreement with the numerical simulations carried 
out. 
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I. INTRODUCTION 

Mathematically, optical dark spatial solitons (DSSs) 
are exact solutions of the one-dimensional nonlinear 
Schrodinger equation (NLSE) for negative nonlinearity 
and nonvanishing boundary conditions 0. Physically, 
they form on background beams of finite width as self- 
supported intensity dips due to the counterbalance be- 
tween beam self-defocusing and diffraction. The required 
negative nonlinearity of the medium causes an inevitable 
reduction of the beam intensity. Losses, saturation and 
high transverse dimensionality result in nonintegrable 
model equations. Despite certain adiabatic relaxation 
characteristics, the solitary solutions of these equations 
have a large number of characteristics J2| in common with 
the soliton solution of the one-dimensional NLSE and are 
widely denoted by the term "dark soliton" . DSSs are gen- 
erated as dark stripes , whereas the only known truly 
two-dimensional optical DSS is the optical vortex soliton 
(OVS) Characteristic of the phase portraits of these 
beams are a ID 7r-phase jump and an on-axis 2tt helical 
phase ramp, also denoted as edge-phase and screw-phase 
dislocations, respectively. 

Optical ring dark solitary waves (RDSWs) were first 
introduced by Kivshar and Yang M. Their evolution 
was studied in the frame of the adiabatic approxima- 
tion of perturbation theory. The quasi-one-dimensional 
treatment allowed the authors to obtain an expression 
for the RDSWs transverse velocity dR/dz as a func- 
tion of the ring radius and the contrast. In addition, 
a linear stability analysis was performed. The authors 
showed that in the small-amplitude limit these ring dark 
waves are described by the cylindrical Korteweg- de 
Vries equation, which has ring soliton solutions Re- 
cently, Frantzeskakis and Malomed derived new equa- 
tions of evolution for small-amplitude solitary waves on 



a finite background [[7J by means of a multiscale expan- 
sion method applied to the generalized cylindrical NLSE. 
Their long-wave solitary solutions propagate "on top" of 
the continuous-wave solution (i.e. on the background) as 
a dark perturbation in Kerr defocusing media. Depend- 
ing on the background beam intensity, they can appear 
both dark and anti-dark for Kerr-like saturable defocus- 
ing nonlinearity. 

The RDSWs may turn out to be of practical interest 
because of their ability to induce waveguides in which 
multiple signal beams could be guided parallel ||. This 
requires detailed analysis of the dynamics of these soli- 
tary structures. The first experimental generation of 
optical RDSWs ^ was conducted with pure amplitude 
modulation in front of the nonlinear medium (reflective 
micrometer-sized dots illuminated by the background 
beam). Under such conditions, pairs of phases oppositely 
shifted by about tt/2 were measured in each diametrical 
slice of the RDSWs jl0|]. In a subsequent experiment 
RDSWs were generated from odd initial conditions 
by using binary computer-generated holograms (CGHs) 
p"2| . Inside a circle the phase was shifted by tt vs. the 
outlying area by shifting the photolithographically pro- 
duced interference lines by half a grating period. In both 
experiments the phase profiles were measured by using 
the multiple-frame interferometric technique for optical 
wavefront reconstruction jL3|. It was found that along 
the NLM the transverse velocity of the RDSWs born 
from even initial conditions is higher than that of the odd 
ones. Numerical simulations have shown that, besides 
the ring radius the transverse velocity of the RDSW 
can be effectively controlled by initial phase modulation 
of the background beam (inside, outside or on both sides 
of the phase dislocation) or by nonlinear interaction with 
a second coaxial dark formation (RDSW or OVS) Jl4| |. 
In saturable media, the RDSWs transverse velocity is 



lower than that in Kerr NLM. At high saturation levels, 
in formal contradiction to but in actual agreement with 
the one-dimensional case, decay of the dark rings into a 
loop of optical vortices of alternate helicities should be 
expected 

We have measured experimentally the dynamics of ring 
dark solitary waves generated by different methods and 
compared those results with the predictions of the ana- 
lytical method and numerical data. In particular, we an- 
alyzed the non-monotonic transverse velocity of RDSWs 
vs. background beam intensity in saturable NLM. The 
RDSWs are generated in the +1 and — 1 diffraction order 
of the background beams. Comparative numerical sim- 
ulations that account for the estimated contrast of the 
RDSWs and medium saturation are carried out. An in- 
tentional reduction of the magnitude of the phase shift 
encoded in the CGHs was found to have a substantial 
influence on the dark ring radii and the contrast at the 
entrance of the NLM. The differences in the nonlinear 
evolution of the RDSWs enable us to verify qualitatively 
the analytical result of Kivshar and Yang (see Eqs. 6 and 
7 in Ref. 10). 



II. RDSW TRANSVERSE VELOCITY VS. 
INTENSITY 

The experimental arrangement is shown in Fig. la. In 
order to obtain a dark ring of a desired radius and a 
phase jump of a certain magnitude, a single-line Ar + - 
laser (A — 488nm) is used to reconstruct the respec- 
tive CGH. The latter are produced photolithographically 
with a grating period of 18fim and ensure a diffraction 
efficiency of nearly 10% in first order. The nonlinear 
medium is ethylene glycol dyed with DODCI (Lamb- 
dachrome). The background beam, diffracted in first or- 
der with the dark ring nested in, is transmitted through 
a slit placed about 15 cm behind the CGH and is gen- 
tly focused on the entrance of the NLM. After passing 
the desired nonlinear propagation path length (0.5cm to 
8.5cm) the dark beam is partially reflected by a prism 
immersed in the liquid and is projected directly on a 
co-moving CCD array with a resolution of 13/im. Alter- 
natively, the dark beam can be recorded at the exit of 
the NLM. For further evaluation the images are stored 
by an image- capturing card. 

The data presented in this section refer to an absorp- 
tion of a = 0.107cm" 1 at A = 488 wn. In a calibration 
measurement we generated a ID DSS by using a CGH 
of the respective type. The soliton constant la 2 (i.e. 
the product of the background-beam intensity I and the 
square of the dark beam width a measured at the 1/e 2 - 
level) is found to reach its asymptotically constant value 
for input powers P^f « 90mW (Fig. lb, dashed curve). 
Since the saturation of the (thermal) nonlinearity is ex- 
pected to have important consequences on the RDSWs 
nonlinear dynamics |l6| ] , we realized a background-beam 



self-bending scheme similar to that used in M| . The re- 
quired asymmetry is introduced by tilting the prism im- 
mersed in the NLM. This results in different propagation 
path lengths for different parts of the beam. The strength 
of the self-bending effect (Fig. lb, solid curve) measured 
in the near field gives an estimated saturation power 
Psat ~ 35mW. The choice of a saturation model for 
absorptive nonlocal nonlinearity is not trivial [jL8[ , Here 
we use a sigmoidal fit of the type /(/) ~ 7/(1 + I/I sa t) 3 - 
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FIG. 1. a) Experimental setup: (CGH - com- 
puter-generated hologram, S - slit, A - aperture, LI and 
L2 - AR-coated lenses (/ = 80mm), NLM - nonlinear 
medium (ethylene glycol dyed with DODCI), Fl, F2 - fil- 
ter sets, CCD - charge-coupled device camera with 13/xm 
resolution, ICC/PC - personal computer equipped with an 
image-capturing card), b) Power dependences of the back- 
ground-beam self-deflection (triangles) and of the quantity 
la 2 (squares) indicating P sa t » 35mW and P^f » 90mW. 
Solid (dashed) curve - sigmoidal fits. 

When a RDSW is generated by pure amplitude mod- 
ulation of the background beam at the entrance of the 
NLM, its transverse velocity increases monotonically 
with the intensity (up to 2.5P s 1 f ; see Fig. 7 in Ref. @). 
By using the CGH of a black ring (with a phase jump of 
7r), we observe a non-monotonic change of the ring radius 
vs. background beam power for a nonlinear propagation 
length of z = 8.5cm (Fig. 2a). The corresponding data 
for the contrast of the RDSW A 2 = (I — I m in)/Io show 
reversed non-monotonic behaviour (Fig. 2b, solid curve) . 
Generally, this is consistent with physical intuition and 
the theory jl6| that the smaller the ring, the darker the 
solitary wave and the higher the total phase shift. As will 
be discussed later (see Fig. 4, the region on the left-hand 
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side of the vertical dashed line), the dark rings initially 
generated with a 7r-phase shift ("black") become wider 
by diffraction and thus arrive as "dark-grey" at the en- 
trance of the NLM. In our case the contrast is lowered to 
A 2 = 0.88 (see Fig. 2b). We measured the change of the 
RDSW's radius at z = 8.5cm as a function of the laser 
power. An evident minimum is observed for a power 
slightly below P}£ ■ For this value the ring contrast is 
highest. Further insight can be gained by evaluating the 
analytical results for the phase of the solitary wave in the 
models of competing cubic-quintic and threshold nonlin- 
earities (see Eqs. 6-9 in Rcf. ]16|] ). With the data for 
A 2 , both models show well-pronounced minima in the 
RDSW transverse velocity for powers slightly lower than 
Pgof (Fig. 2, dashed curve). It should be noted that the 
experimental data in Fig. 2 were obtained by compen- 
sating losses during the propagation with an appropriate 
change of the input power. However, even for fixed input 
power the RDSW changes its transverse velocity along 
the medium. In that sense Fig. 2a presents the effect of 
the enhancement of the phase shift and the contrast of 
the RDSW by saturating the nonlincarity in an integral 
form. This result, however, does not explain why in our 
previous experiments only a monotonic increase of the 
RDSW radii was observed |9) . 
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FIG. 2. a) Power dependence of the RDSW radius after a 
nonlinear propagation path length z = 8.5cm. b) Power de- 
pendence of the RDSW contrast (dots) and estimated trans- 
verse velocity (dashed curve) for saturable cubic-quintic non- 
linearity. See text for details. 

In order to explain this discrepancy we look at the 



evolution of the RDSW before it enters the NLM. Nu- 
merical simulations are carried out for the evolution of 
"grey" solitary waves with initial radii Rq = R(z = 0) 
twice as large as the ring width tq = r[z = 0). The nu- 
merical procedure used to solve the 2D NLSE is based on 
the beam propagation method over a 1024 x 1024 grid. 
In agreement with our measurement, initially "black" 
rings generated by CGHs become "grey" after passing 
four Rayleigh diffraction lengths Lmff = kr^ and reach 
the estimated contrast A 2 = 0.88. The evolution of the 
RDSWs radius for different intensities is calculated for 
up to 10 nonlinear lengths (Latl = (fc|^2|^o) _1 , k is the 
wavenumber, and |ri2|Jo the nonlinear refractive index 
correction). The behaviour of the RDSW radius in the 
NLM can be summarized as follows: 

i) For a fixed input intensity in the range of (0.2 — 2.2)1^ 
all RDSWs have their minimal ring radius for propaga- 
tion distances between 1.5-Ljvl and 3Ljvl- 

ii) For nonlinear propagation distances shorter than 
2L jvl the grey rings monotonically reduce their radii for 
increasing intensity (Fig. 3, lowest box). 

iii) For long propagation distances (> 6Ljvl) this ten- 
dency is reversed (Fig. 3, upper box). 

iv) In between, the radius of the RDSW as a function 
of the intensity shows a characteristic minimum around 
4.5Ljvl (Fig. 3, middle box). This propagation distance 
is in reasonable agreement with the NLM length of 8.5cm 
corresponding to ALnl, for which Fig. 2 is obtained. 
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FIG. 3. RDSW radius vs. normalized background-beam 
intensity at three characteristic nonlinear propagation dis- 
tances. Model parameters: Ro/ro = 2, A 2 = 0.88, s — 0.4. 

Generally, we found qualitative confirmation of the ex- 
istence of a minimal radius R m in of the RDSW, which 
depends on the initial contrast A 2 [| 16 1. The evolu- 
tion of the contrast itself is strongly influenced by the 
saturation of the nonlinearity. In view of the above it is 
interesting to recall the observation which we made when 
modeling the decay of multiple-charged OVSs under ra- 
dial perturbations [fl9f : The small-amplitude even-phase 
initial conditions for perturbation and the saturation of 
the nonlincarity resulted in a negative transverse velocity 
of the rings. Their contrast was found to increase with 
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decreasing radius and the radial phase variation reached 
a maximum of 0.757T at 5Ltvl- The initially decaying and 
still strongly overlapping optical vortices partially recov- 
ered their multiple-charged state (Fig. 7 in Ref. fl9| ). Pe- 
riodic sign changes in the transverse velocity of the dark 
rings in the propagation direction were observed in the 
simulations up to IOLnl with a period of approximately 
2Latl- 



III. RDSW TRANSVERSE DYNAMICS VS. 
INITIAL PHASE SHIFT 

It should be pointed out that initially "black" RDSWs 
with high transverse dynamics are not likely to be present 
in a real experiment, at least when generated by holo- 
grams. The reason lies in the necessity to select one of 
the diffracted beams, which requires a certain free-space 
propagation. After some diffraction is "accumulated" , 
the dark rings with small radii spread out substantially 
and their contrast decreases. Dark rings with large radii 
are less affected and they evolve more slowly inside the 
NLM. Figure 4 is intended to visualize this before and 
inside the NLM. The coordinate z = is fixed on the 
nonlinear interface. Solid and dashed lines correspond to 
Kerr and saturable Kerr nonlincarity (s = I /I sa t = 0.4). 
The simulation presented is carried out for / = 7 S( ^ and 
for an initial phase jump of ir, which will be flattened 
in the course of propagation. Generally, the saturation 
of the nonlincarity leads to a reduction of the dynamics 
of the RDSW. For our further considerations it is im- 
portant to note that there is no difference whether the 
phase inside the ring is bigger or lower than the phase 
<& „t outside the ring, provided A<& = |$ in — <& ont | = 7r. 
Equivalently, there should be no difference in the trans- 
verse dynamics of the RDSWs generated by the +lst and 
— 1st diffraction order beam when a phase jump of n is 
encoded in the CGHs. This was carefully proved by nu- 
merical simulations and experimentally. 
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FIG. 4. Diffraction-governed evolution of odd RDSWs with 
initial 7r-phase jumps over 4 Rayleigh diffraction lengths 
followed by nonlinear propagation in a Kerr (solid curves) 
and saturable Kerr-like medium (dashed curves; s = 0.4). 



A. Linear evolution for |A$| < tt 

The situation changes when we produce dark ring 
waves encoded on CGHs with phase shifts smaller than 
it. Unfortunately, in binary CGHs only ir/N steps are 
possible, where one grating period consists of 27V elemen- 
tary stripes. In our case N = 3 and the grey rings are 
designed for |A$| = 2ir/3. The beam of a single-mode 
He-Ne laser is expanded 3.5 times and reproduces the re- 
spective CGH. By means of a lens (/ = 4cm) the beam 
transmitted directly and the ±lst order diffracted beams 
are imaged on a screen located about 100cm behind. By 
varying the lens position we are able to follow the dark 
ring evolution. Figure 5a shows the images of the still 
overlapping first-order diffracted beams, which interfere 
with the zeroth order one. It can be seen clearly that 
the rings with |A$| = 2tt/3 have smaller (larger) radius 
than each of the dark rings with an initial phase shift 

|A$| = 7T. 



a) 



b) 32 





li 



AO=2ti/3 



AO=7T. 




I = It 



Vertical dashed line - the nonlinear interface. 



FIG. 5. Linear propagation of ring dark beams: a) In- 
terference patterns indicating different transverse velocities 
of free-propagating, initially grey (|A$| = 2-7r/3) and black 
(|A$| = 7r) ring dark beams at L = 5.5cm. b) Host graph: 
ring radii vs. CGH-to-CCD-array distance L. Inset: the 
same vs. lens-to-CGH-distance. Dots - List diffraction order 
beams, blank circles - initially black dark rings. See text for 
details. 

The ring radii observed for lens positions from 4cm to 
16cm are shown in the inset of Fig. 5b. In order to get 
better spatial resolution, the measurement was repeated 
with a CCD array illuminated directly after the CGH 
(Fig. 5b). Blank circles denote the dark ring radii for 
|A$| = 7r, the solid circles refer to the ring with encoded 
phase jump | A$| = 27r/3. Note that in the two graphs of 
Fig. 5b the ring radius R is denoted in CCD-camera pix- 
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els, but these units are not directly comparable because 
of the different focusing conditions. On the CGHs the 
ring radius is encoded in 20 pix. (i.e. Rq = 60/xm). It is 
easy to understand the observed differences: The diffrac- 
tion tends to broaden the dark rings and flatten the phase 
jumps. As a result, for <F;„ < $ OM t an effective concave 
phasefront is formed which "focuses" the dark ring. In 
the other diffraction order (fr ou t appears to be bigger than 
$i n . The phase change induced by the diffraction results 
in an effective convex phasefront that increases the dark 
ring radius. For |A$| = 71", the diffraction influences both 
first-order beams in the same way, but more weakly. 



B. Nonlinear evolution of initially grey RDSWs 

Although the nonlinear regime was analytically ana- 
lyzed by Kivshar and Yang || , to the best of our knowl- 
edge it has not been investigated experimentally. The 
basic results (Eqs. 5-7 in Ref |(|) can be summarized as 
follows: 

i) For each RDSW there is a minimum ring radius R m in 
for which the RDSW has its highest contrast. 

ii) Depending on the initial value of the phase jump, the 
RDSW either collapses along the NLM to reach R m in and 
diverge monotonically thereafter, or it diverges monoton- 
ically from the beginning. 




0123456789 10 
z/L NL 

FIG. 6. Nonlinear evolution of initially grey RDSWs in a 
Kerr (solid curves) and saturable Kerr-like medium (dashed; 
s = 0.4) for positive (a) and negative (b) phase shift 
A$ = $ 4n — $ ut- The model parameters correspond to 
A 1 = 0.88 (i-e. to |A$| = 0.8377T) and / = 1^ . 



Figure 6 is intended to clarify these conclusions. We 
solved the NLSE numerically by keeping the following 
parameters close to the experimental values: RDSW 
radii ((R/r)\ z= o = 2,4, and 6), background-beam in- 
tensity (/ = IgoY), initial contrast A 2 = 0.88, phase 
jump |A$| — 0.837-7T, and saturation of the nonlincar- 
ity s — 0.4. In Fig. 6a the corresponding results for the 
case of a higher phase inside the rings (<!>„ > § ou t) are 
presented. In this case, the RDSWs start diverging im- 
mediately after entering the NLM. The data shown in 
Fig. 6b refer to the case of a lower phase inside the dark 
ring and higher phase outside (<&i ra < $ ou t)- The ef- 
fective concave wavefront forces the RDSWs to collapse 
to a minimum ring diameter. Thereafter they start to 
broaden faster. The main tendency, that the broader the 
RDSW initially is, the longer the nonlinear propagation 
needed to reach Rmin will be, is well expressed. The 
other general tendency, that the smaller the RDSW di- 
ameter is, the higher its transverse dynamics is, can also 
be clearly seen. 

Both tendencies and the influence of the saturation 
seem to be understood theoretically, but still need ex- 
perimental confirmation. We performed additional mea- 
surements with the experimental setup shown in Fig. 1, 
in which the telescope is replaced by a single lens 
(/ = 12cm) in order to focus the beams near the en- 
trance of the NLM. Moderate and strong saturation s 
can be achieved by raising the absorption a to 0.2cm -1 
(P^f = 25mW;P sat = 40mW; s = 0.6) and 0.4cm" 1 
(P^f = 20mW; P sat = U{±5)mW; s w 1.4). In order to 
keep the conditions for both first-order beams as similar 
as possible, the beam transmitted directly was blocked 
at the entrance of the NLM and both ±fst order beams 
were recorded simultaneously. Unfortunately, due to the 
interaction of the two background beams, the interaction- 
free nonlinear propagation path length is limited. The 
results of the measurements are shown in Fig. 7a for 
moderate and high saturation (host graph and inset, re- 
spectively). Blank circles denote "black" RDSWs radii, 
solid ones RDSW radii for A$ = ±27r/3. The best re- 
sults for moderate saturation are achieved with CGHs 
containing phase dislocations with radii of R = 30/im. 
For high saturation, rings with dislocation radii of 60/im 
are used. As seen, after some 5mm of nonlinear propaga- 
tion the "grey" RDSWs reach R m i n and diverge mono- 
tonically afterwards. The beam of the other diffraction 
order and the RDSW generated for |A$| = it start di- 
verging from the beginning. The interpolated curves in 
the figure are intended to guide the eye. The interpo- 
lated curves for z < 2mm and z > 17mm are not shown 
for an obvious physical reason: As mentioned, the grey 
RDSWs do not enter the NLM with equal radii and equal 
(absolute) phase shifts. For larger propagation distances, 
the saturation-like behaviour is caused mainly by interac- 
tion with the second co-propagating beam. The reduced 
transverse velocities of the RDSWs for higher saturation 
are also strongly pronounced. The grey-scale images in 
Fig. 7b show the RDSWs at z = 2mm and 8mm inside 



5 



the NLM. 
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FIG. 7. Nonlinear propagation of RDSWs: a) Evolution 
of the RDSW radii along the nonlinear medium at moderate 
(host graph, s = 0.6) and high saturation (inset, s ~ 1.4). 
Solid and blank circles - initially grey (|A$| = 2n/3) and 
black (|A$| = 7r) RDSWs, respectively, b) Grey-scale images 
of RDSWs at z = 2mm (upper row) and 8mm (lower row) in 
a moderately saturated NLM. 




FIG 8. Gray-scale images showing the development of an 
internal substructure of a RDSW (Ro = 90/j.m) at moderate 
saturation. Left to right: z — 2,5, 8, 11, 14, and 17mm. 



posite quasi-2D phase dislocations. The relation of these 
coaxial rings to the small-amplitude dark solitary waves 
described by Frantzeskakis and Malomed should be 
clarified. 



IV. CONCLUSION 

It is worth noting that it is difficult (if possible at all) to 
ensure perfect odd initial conditions (radial phase jumps 
equal to 7r and contrast equal to unity) for the genera- 
tion of RDSWs in cubic nonlinear media. When such a 
phase jump is encoded in a CGH, there is no difference 
between the evolution of the two beams in first diffraction 
order. This changes gradually for smaller phase jumps. 
In qualitative agreement with the theory Q we presented 
evidence of the monotonic increase of the dark ring radius 
for > $ ou t and the corresponding non-monotonicity 
in the opposite case. For grey initial conditions and sat- 
uration of the nonlinearity we found characteristic re- 
gions for the dependence of the RDSW radius on the 
background-beam intensity. The radius can increase, de- 
crease monotonically or go through a minimum. The 
results presented are in qualitative agreement with the 
analytical ones regarding the phase of the solitary wave 
in the presence of saturation Pq| . 
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By evaluating the data corresponding to the largest 
RDSWs in our experiment (R — 90/im) we observed 
internal ring splitting for moderate saturation (Fig. 8). 
Despite their large radii, these RDSWs have non-zero 
transverse velocity. They broaden and flatten due to the 
saturation (see Fig. 8, frame at z = 2mm) while evolving 
into grey rings. It seems that such a grey ring can serve 
as a background on which two coaxial rings can form, 
subsequently narrow and thus get a high modulation 
depth for increased propagation distance (Fig. 8, frames 
at z = 5, 8, 11, 14, and 17mm). The mechanism for their 
formation includes effects from increasing a phase shift 
and decreasing a transverse velocity due to the saturation 
of the nonlinearity J16| , as well as repulsion between op- 
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